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Abstract 

We study a two-dimensional gas of inelastic smooth hard dimers. Since the collisions between 
dimers are dissipative, being characterized by a coefficient of restitution a < 1, and no external 
driving force is present, the energy of the system decreases in time and no stationary state is 
achieved. However, the resulting non equilibrium state of the system displays several interesting 
properties in close analogy with systems of inelastic hard spheres, whose relaxational dynamics has 
been thoroughly explored. 

We generalise to inelastic systems a recently method introduced [G.Ciccotti and G.Kalibaeva, 
J. Stat. Phys. 115, 701 (2004)] to study the dynamics of rigid elastic bodies made up of different 
spheres hold together by rigid bonds. Each dimer consists of two hard disks of diameter d, whose 
centers are separated by a fixed distance a. By describing the rigid bonds by means of holonomic 
constraints and deriving the appropriate collision rules between dimers, we reduce the dynamics 
to a set of equations which can be solved by means of event driven simulation. After deriving the 
algorithm we study the decay of the total kinetic energy, and of the ratio between the rotational 
and the translational kinetic energy of inelastic dimers. We show numerically that the celebrated 
Haff's homogeneous cooling law t" 2 , describing how the kinetic energy of an inelastic hard sphere 
system with constant coefficient of restitution decreases in time, holds even in the case of these non 
spherical particles. We fully characterize this homogeneous decay process in terms of appropriate 
decay constants and confirm numerically the scaling behavior of the velocity distributions. 

PACS numbers: 02.50.Ey, 05.20.Dd, 81.05.Rm 
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I. INTRODUCTION 



Over the past two decades much progress has been made in understanding the dynamics 
of granular matter, i.e a collection of macroscopic particles interacting via a short range 
repulsive potential in which energy is lost in inelastic collisions^**- The industrial and 
technological importance of products which are either powders or granulates is now becoming 
generally recognized. Such a continued interest has unveiled a series of surprising phenomena 
and determined new challenges to physicists, because of the simultaneous presence of many 
body effects and non thermal fluctuations. In fact, in granular materials such as sand, cereals, 
steel balls etc. the ordinary temperature is irrelevant and consequently thermal equilibrium 
is not realized. Due to inelastic collisions a granular system, when not driven cools, i.e. 
the one particle velocity distribution tends to narrow in time 4,5 '. The inelastic hard sphere 
(IHS) model represents, perhaps, the simplest description of a granular system and most 
of the theoretical investigations are based on it. A collection of IHS sufficiently rarefied 
is called granular gas, because its dynamics consists of ballistic motions alternated with 
instantaneous binary collisions which conserve total momentum. The inelastic collisions are 
modeled by means of constant coefficient of normal restitution a, which determines the rate 
at which the kinetic energy is dissipated. Therefore, we possess nowadays a large amount of 
information concerning the IHS which has become the standard reference model for granular 
matter e^WV?, 

A natural question to ask is whether the behaviors observed in the IHS are generic or 
peculiar. The relevance of such an issue is not academic, since many actual materials, such 
as rice, needles, are far from being spherical but are elongated. As far as thermal systems are 
concerned the role of the shape of the constituent molecules in determining their thermody- 
namic behavior has been recognized since the pioneering work of Onsager™. Entropic effects 
(packing) of sufficiently anisotropic molecules (prolate ellipsoids, needles, sphero cylinders, 
disks) determine the formation of mesophases 1 ^. In granular matter, instead, the effect of 
the non spherical shape of the particles is still largely unexplored in spite of its potential 
relevance in many phenomena such as nematic ordering transitions or violations of the en- 
ergy equipartition. Few exceptions to such a trend are represented by the investigations of 

rc f„ 15.16.17.18 

Since the computer simulation is a fundamental tool in the study of granular systems, it 
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is crucial to devise an efficient method to treat the collisions, which usually represents the 
most time consuming part of a numerical code. The Event Driven (ED) method reduces such 
an effort to the calculation of the collision times, since the interactions in hard core systems 
occur only at the moment of collision. Once the shortest collision time has been computed 
all the particle positions are propagated inertially with the new velocities determined by 
the collision rules. The advantage over the fixed time step method in the case of hard 
objects is clear since there is no algorithmic error in the integration of the equations of 
motion. However, in the case of non spherical bodies the calculation of the collision time 
is not simple. Therefore, some authors resorted to a fixed time step method. Recently, 
Ciccotti and Kalibaeva^ have shown that it is possible to include a constraint, such as 
a fixed bond length between two elastic hard spheres, and preserve the advantages of the 
ED simulation. Hereafter we illustrate how to extend the ED method to a model of bond 
constrained inelastic hard disks. 

The structure of the paper is the following. In section II we define the model and derive 
the equations of motion for a pure system of inelastic dimers. We also generalise the method 
to mixtures of inelastic dimers and disks and to systems of dimers in the presence of a fixed 
impenetrable wall. In section III we present numerical results obtained by means of the ED 
simulation. We discuss the properties of the homogeneous cooling state of inelastic hard 
dimers and obtain numerically the relevant parameters and the distribution functions. In 
section IV we draw the conclusions and discuss the future developments. 

II. MODEL AND ALGORITHMS 

In the present paper we discuss the properties of a system comprised of N dimers con- 
sisting of two identical hard disks of mass m rigidly connected. The diameter of the disks 
is d and the distance between their centers is a. In Fig. ^ there is a snapshot of the system. 
By describing the rigid bond between the two disks via an holonomic constraint we are 
able to apply the ED numerical simulation technique. 

Let us consider the two-dimensional motion of an assembly of rigid dimers. We make the 
assumption of rigid body dynamics: 

• The duration of contact is negligible and the interaction forces are high, so the velocity 
changes are nearly instantaneous without notable change in positions. 
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• The contact area is also negligible and the deformations are small in the impact zone 
so that the impact occurs at a single point of each body. 

• Double contacts cannot occur simultaneously. 

The motion of each dimer, A, is described by the velocity of its center of mass and by 
its angular velocity uj a . The discontinuities of linear and angular momentum occurring at 
a collision event are obtained by imposing the conservation of the total linear and angular 
momentum and the law governing energy loss. 



A. Free streaming 

In between collisions the dimers in the absence of external forces perform an unperturbed 
roto-translational motion. As shown in ref.— such a free streaming can be described by the 
following parametric equations for the coordinates rj.(t) and r 2 (£) of the centers of the two 
disks forming the dimer: 




FIG. 1: Snapshot of a system of 400 dimers with a = d and a = 0.95. 



v x (t) = R A (0) + V A (0)t - - |r 2 (0) - n(0)J cos(u A t) - 
r 2 (t) = R A (0) + V A (0)t + - [r 2 (0) - n(0)] cos(aMt) + 



v 2 (0)-vi(0) 

2u A 
v 2 (0) -vi(0) 

2u A 



sin(uj A t) 



sin(ujAt) 



(1) 
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where Ra(0) e V^(0) are respectively the position and the velocity of the center of mass 
of the dimer A at the instant t = 0. The first two terms describe the translational motion, 
while the two last describe the rotational motion about the center of mass of the dimer. The 
angular velocity uja is obtained by solving the dynamics in the presence of the holonomic 
constraint 

X =[r 2 (t)-n(t)] 2 -a 2 = (2) 
In fact differentiating with respect to t eq. (J2J) we obtain: 

! cos(ujAt) | 



[v 2 (t) - v x (*)] ■ [r 2 (t) - ri(t)] = -u A \ [r 2 (0) - n(0)] sm{u A t) - [v 2 (0) - Vi(0)]- 



•|[r 2 (0) - n(0)] cos(w A t) + [v 2 (0) - Vl (0)] Sm(a;At) | = 

(3) 

Using eqs. (J2J) and (|3J) we obtain that angular velocity is constant and given by: 



B. Collision times 



In order to give a self-contained account of the method we include the determination of 
the collision times following the presentation given in ref.— . The condition for the collision 
between two "molecules" A and B, formed by the four atoms 1, 2 and 3, 4, respectively are: 





-ri(t)} 2 


= d 2 


Ht) 


-r 2 (t)] 2 


= d 2 


Ht) 


-n(t)} 2 


= d 2 


Ht) 


-r 2 (t)] 2 


= d 2 



(5) 



Substituting the equations (0) for the position of each hard disk as a function of time in the 
collision conditions (0) we get a set of equations in t, from which we take the smallest one 
as the collision time of the two molecules. 
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B.l Dynamics of Collisions 



In each collision between two dimers four disks are involved, therefore we write the 



conservation laws for the total linear momentum 

m(v lf + v 2/ + v 3/ + v 4/ ) = m(v H + v 2i + v 3i + v 4i 



(6) 



and for the total angular momentum 

m(ri x Vi/ + r 2 x v 2/ + r 3 x v 3/ + r 4 x v 4/ ) = m{r x x v H + r 2 x v 2i + r 3 x v 3i + r 4 x v 4i ) (7) 

where the subscripts i and / indicate pre-collisional and post-collisional observable, respec- 
tively. Since the variables r 3 - (j = 1,4) do not change during the instantaneous collision, we 
do not need to distinguish between their pre-collisional and post-collisional value. 

Upon differentiating the constraint (J2J) we obtain relations between the velocities of the 
particles belonging to the same dimer: 



(V 2 j 


- v u ) 


ri2 


= 


(v 2/ 


-Vi/) 


ri2 


= 




- v 3i ) 


r 3 4 


= 


k ( V 4/ 


- v 3/ ) 


r 3 4 


= 



(8) 



with Tjk = Tj — r^. We simply satisfy eq. (jUJ) by writing the change of momenta of the two 
molecules A and B with the help of a vector Av, to be specified in the following 



m(vi/ + v 2/ ) = m(vii + v 2i ) + mAv 
m(v 3f + v 4/ ) = m(v 3i + v 4 j) - mAv 



(9) 



In order to satisfy simultaneously eq. (JUJ) and the constraints (jSJ) it has been shown, assuming 
that the colliding particles are 1 and 3, that one obtains the following expressions^ 



Vi/ = Vu + Av - 



Av ■ r 



12 



2a 2 



-ri2 



Av • r 



v 2 f = v 2i 
v 3/ = v 3i 
v 4/ = v 4i 



12 



2a 2 



-ri2 



Av H — — r 34 



(10) 



Av ■ r 34 

2a 2 



2a 2 



-r 34 
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The change of the angular momenta of the two colliding dimers are 

m(r 1 x vi/ + r 2 x v 2/ ) = m(r 1 x vu + r 2 x v 2i ) + rari x Av 



(11) 



™(r 3 x v 3/ + r 4 x v 4/ ) = m(r 3 x v 3i + r 4 x v 4i ) - mr 3 x Av 

so that the conservation of the total angular momentum at collision implies that Av is 
directed along the direction of the vector ri 3 connecting the centers of the two colliding 



disks, i.e 



Av = Av 



ri3 



(12) 



Finally substituting eq. (fT2"j) into eq. (fTUjl we arrive at 



r 43 i"i 3 • i°i 2 
Vij + Av— Av r i2 

d 2a 2 d 



a r 13 • r 12 

v 2 / = v 2i + Av - - - r 12 



r 3f 



2a 2 d 
r 43 rx 3 • r 34 
v 3i — Av— - + Av 7r-—r-iA 

3 d 2a?d 6 



(13) 



i"i 3 • r 34 

v 4/ = v 4i - Av r 34 

In order to fix the amplitude Av we consider the change of total kinetic energy occurring in 
a single collision event: 

T 2 , r 2 , r 2 > ,f ' , , r 2 ,,.2 ,,2 , r 2 . 

(14) 



AE = j( V h + V 2/ + V 3/ + Vl/) - y K + V 2 2 , + V 3 2 4 + V 2 i; 



" I^d? [(ri3 ' 1-12)2 + (ri3 ' r34)2] } 



mAv 
d 



(vsi - v u ) • r 
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For a perfectly elastic collision the condition AE = determines Av: 



Av 



(vh - v u ) • r 13 /rf 



(15) 



(ri3 ■ r 12 ) 2 + (ri 3 • r 34 ) 2 

To generalize the collision rule to the case of inelastic collisions we consider the variation of 
the relative velocity of the two colliding disks along the center to center direction: 



(v 3/ - Vi/) • n 3 = (v 3 , - vij) • r 13 - 2Afrf|l - j^-p [(rig • r 12 ) 2 + (r 13 • r 34 ) 2 ] J (16) 

and impose the condition that the projection of the relative velocity along r 13 after the 
collision is proportional to the projection of the relative velocity before the collision via a 
coefficient of restitution a (0 < a < 1), which takes into account the energy dissipation. 



(v 3/ - Vif) ■ r 13 = -a(v 3i - vi,) ■ r 
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(17) 



In this case one sees that the rule (|15|) is modified 

= 1 + a (v 3 » - Vh) • ri3 

2rf 1 _ ^ [ (ri3 . ri2 )2 + ( ri3 . r34 )2] ^ ^ 

By substituting the relation (jl8)l into eq. (fPfj) we find the corresponding energy loss in 
an inelastic collision 

AF-- l ~ a2 [(V3» ~ v^) -ri 3 ] 2 

Let us compare the energy dissipated (fT^j) by two inelastic dimers with that associated with 
the collision of two inelastic hard disks of diameter d, mass m and coefficient of restitution 
a: 

I - a 2 9 
AE dd = -™—£$r [(vm - v H ) • r 13 ] (20) 

and notice that in the first case the amount dissipated depends on the relative orientation 
of the two dimers, the largest dissipation occurring for head to head collisions, where the 
two axes v\2 and are parallel to 1*13. This orientational dependence of the dissipation of 
dimers renders the model similar to a model of inelastic disks with a fluctuating coefficient 
of restitution. 



B.2 Collision between a dimer and a disk 

he behavior of mixtures of different grains is also a subject of extreme interest. The 
dimer-disk mixture has not been studied so far. It is a simple exercise to derive the collision 
rules in the case of the collision between a dimer and a hard disk having the same coefficient 
of restitution, 

/ A _ _ 

12 

— ri2 



Vi/ 


= Vh + 


Au ■ r 
Au — 

2a 2 


V 2 f 


= v 2i + 


Au ■ ria 

2a 2 ri2 




= v di - 


m . 
-Au 

M 



(21) 



where we have assumed that the mass of the disk is M. is the position of disk center 
and a its diameter. We also introduced the vector Au: 

a. ' V «7"L (R '7 ) /V < 22 > 

1 + f - 2^K R rf- r l)- r 12j G 
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with C = (d + cr)/2 and is the coefficient of restitution. 

In this case the total energy dissipation in a dimer-disk collision is 

Aft— .^g [CV«-Vu) •(*-*,)]■ (23) 

B.3 Collision between a dimer and a smooth impenetrable wall 

second extension of the previous methods concerns the collision of a dimer against an 
impenetrable hard wall with a restitution coefficient a w . Taking the limits a — > oo, M — > oo 
and = 0, we find in this case: 

Aw ■ r 12 



Vi/ = v u + Aw — — ri 2 

Aw • ri2 

v 2 / = v 2i H — -5 — ri2 

2a 2 



(24) 



with 



Aw = -(l + aM ) ? rfi n fl (25) 



and n is a unit vector normal to the surface of the wall. 
The relative total energy change is 



A ^ " ~ m 2 l-^(n.r 12 ) 2 (26) 



III. COOLING DYNAMICS 



In the present section we focus on the evolution of an initially homogeneous gas of inelastic 
(a < 1) hard dimers constrained to move on a plane. The system consists of N particles 
contained in a square subject to periodic boundary conditions and is not externally driven. 
Inelastic collisions between the dimers result in a loss of kinetic energy and thus determine 
a "cooling" of the gas. 

This problem has been attracting much interest since the inelastic collisions make the 
system to behave very different from standard molecular fluids and plays a central role in 
this sub-area of non equilibrium statistical mechanics. The discovery of the homogeneous 
cooling state (HCS), characterized by a one particle distribution function whose shape is 
self similar, if the velocity is appropriately rescaled, has triggered a lot of attention 20 . 



10 



In order to put things in perspective we take as a reference system, a set of N inelastic 
hard disks of diameter a in two spatial dimensions (D = 2), since its properties have been 
thoroughly studied. Such a system starting from an equilibrium state of a corresponding 
elastic hard disk system, cools down uniformly, and decreases its average kinetic energy per 
particle Ek in (t) = T(t) according to Haff's law^ 

Tit) = - - (27) 

K) (l + 7 ^) 2 1 ' 

where T(t) is the so called granular temperature, T the initial temperature, z/ is the equi- 
librium Enskog collision frequency at T and the non dimensional parameter 7 is: 

Remarkably the homogeneous cooling law holds in any spatial dimension, D, as far as smooth 
inelastic hard spheres are considered. 

HCS is characterized by uniform density and velocity fields and is maintained only for a 
finite time, or for a finite number of collisions, since it is linearly unstable to fluctuations. 
After the initial homogeneous stage, the system enters an inhomogeneous cooling state where 
a vortex structure develops in the velocity field (shearing instability) followed by a density 
instability (clustering)^*. 

A theoretical estimate of z/ for elastic hard disks is provided by the following Enskog 
expression 

v = 2^na^-g(a) (29) 

where n is the number of particles per unit area and g (a) is the hard disk equilibrium pair 
correlation function evaluated at contact: 

9i?) = ( 30 ) 

and T) is the area fraction 77 = 7mer 2 /4. It is also useful, for future comparison, to express 
the number of atomic collisions per particle, r, suffered by the particles until time t by the 
phenomenological relation connecting to the inelasticity parameter 7 and the equilibrium 
collision frequency: 

r(t) = ^^ = -ln(l + 7*bf) (31) 
1\ 7 
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A second important property of the HCS concerns the velocity distribution function, 
which of course cannot be stationary due to the energy loss, but after a short transient 
assumes a form which depends on time only through the granular temperature, T(t) = 

«'■*)" SSF*^) (32) 

with 

vl{t) = ih J dDvv2 K v ^ ^ 

The scaling function, $(z) is time independent in the HCS and in the limit of small 
dissipation approaches a Gaussian, i.e. Q>(z) ~ ix~ D l 2 exp(— z 2 ). For smaller values of a a 
perturbation theory about the elastic state accounts for the departure from the Maxwellianii. 

Is the scenario described above preserved in a gas of hard inelastic dimers? In other 
words can we find constant values and v^f* such that eqs. (|27|) and (|31|) still hold 
with these effective parameters? Secondly, do the distribution functions of the translational 
and rotational components have the scaling property (|32)1? If the answer is affirmative the 
inelastic hard dimer system possesses a homogeneous cooling state, whose importance for 
the study of granular gases has been stressed by various authors 20.21.22.23.24.25 _ 

The evolution of the inelastic dimer system can be characterized in terms of the average 
values per particle of the total kinetic energy, the translational and the rotational kinetic 
energy. The last two quantities averaged over the particles define the partial granular tem- 
peratures: 

N 

T tr (t) =< E tr (t) >= ^ E t V ^) " < V ^)>] 2 (34) 



A=l 



T rat (t) = 2 < E rat (t) >= (35) 

A=\ 

where V^(t) is the center of mass velocity of a single dimer A, uj^ its angular velocity and 
I = ma 2 /2 its moment of inertia. Our results for the decay of the average total kinetic 
energy E(t) = T tr (t) + \T rot {t) are shown Fig. El We observe that the cooling process occurs 
according to the same universal inverse power law t~ 2 which characterizes the inelastic hard 
sphere systems. In addition, when plotted as a function of the total number of atomic 
collisions suffered on the average by a dimer times the inelasticity parameter (1 — a 2 ), the 
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curves corresponding to different values of a nicely fall one onto the other. However, the 
slope characterizing the various dimer systems is different, being in fact slower, from the 
corresponding slope of the hard disk systems. This feature represents the fingerprint of the 
structure of the dimers. 

A second remarkable feature is the validity of the logarithmic law (|31|) (see Fig. |3J), 
with the effective parameters i/g and = (1 — a 2 )/(2D e ) shown in table I. Notice that 
D e = D in the case of inelastic D-dimensional smooth hard spheres. We note that while D e 
is nearly independent of the elongation, a, of the dimer, the initial collision frequency i/q 
does depend on a, being connected with the total cross section of this object. 




FIG. 2: Data collapse for dimers having different coefficients of restitution obtained by plotting the 
decay of total energy as a function of (1 — a 2 )N atco u/N, i.e. the total number of collisions suffered 
by a molecule multiplied by the inelasticity factor (1 — a 2 ). The dashed line corresponds to the 
energy decay of the reference hard disk system and shows a faster decay. In the inset the same 
quantities are shown as functions of the rescaled time VQ^t, showing the expected Haff's inverse 
power law decay. 

The entries of table I have been obtained by fitting the numerical data with the logarith- 
mic law eq. ()31)1 with D e and used as fitting parameters. Our data analysis reveals that 
in all cases considered the fit is very good and that D e is very close to 3. 

We use a simple heuristic argument in order to derive the slower cooling rate of the 
hard dimers with respect to the disks. Let us compare formulae (|19|) and ()20|) . representing 
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FIG. 3: Average number of atomic collisions per dimer rescaled by the inelasticity factor (1 — a 2 ) as 
a function of the rescaled time for various values of the inelasticity and for a center-to-center 
distance a = d. 

the energy dissipated in a collision between two dimers and two disks, respectively. What 
renders the behavior of the two systems different, when the coefficient of restitution a of 
both is the same, is the presence of the angular weighting factor, /: 

1 



./ 



2-i[(r 13 T 12 )2 + (r 1 3T 3 4) 2 ] 



(36) 



We can rewrite eq. (fTHjl as 



I- a 2 

AE = -m 2rf2 / [(vai - vii) ■ r 13 ] ' 



(37) 



We consider the average value of eq. (J37j) over many collisions, here indicated with the 
symbol < ■ > and assume the following factorization 

1 - a 2 , 1 



(AE) = -m^— (/>-=( [(v» - v u ) ■ r 13 n 



2 w ' d 2 

The last average turns out to be proportional to the average total energy, i.e.: 



(38) 



1 ,r, l2 . (E) 



m 



which gives 



(AE) = - 1 -^(f)(E). 



(39) 



(40) 



14 



We compute (/), assuming a dimer elongation a = d. It is straightforward to see that the 
angles formed by the directions ri3 and r±2 or are subjected to the constraint, required 
by the condition of non overlapping two dimers, that only the angles comprised between 
— 2tt/3 and 27r/3 are allowed (see Fig. 0J). We further assume a uniform distribution of 
these angles in the above interval, a guess which is confirmed numerically. As shown in Fig. 
the computed distribution is almost uniform. Using this fact, we perform the average 
of formula (j3T?j) obtaining the value (/) = 0.64. With this value we obtain the effective 
"dimension" D e = D/(f) = 3.12 which is compatible with the numerical values shown 
in table I. Intuitively this result is consistent with the idea that the energy stored in the 
rotational degree of freedom is only indirectly affected in the dissipation process and provides 
a kind of reservoir slowly dissipated. 




FIG. 4: Allowed angle between the directions ri3 and r±2 when two dimers are at contact. Since 
the dimers cannot overlap only the angles between — 2tt/3 and 2ir/3 are allowed. 



A. Energy equipartition 

One of the most peculiar features of granular systems is the lack of energy equipartition 
among different kinetic degrees of freedom. Such a phenomenon has been observed both in 
driven and undriven systems of smooth and rough hard-sphere mixture o 26 ' 27 ! 28 ' 29 and during 
the homogeneous cooling of a gas of inelastic needles^. 
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FIG. 5: Distribution of the angle between ri2 and ri3 measured in subsequent intervals, each 
containing 10 4 collisions. The distribution has been obtained for 900 dimers, a = 0.98. The elastic 
case a = 1.0 is also shown (bullets) for comparison. 

In the case of dimers, the translational and rotational degrees of freedom satisfy the 
equipartition of energy only when a = 1, as required by equilibrium statistical mechanics. 
We have evaluated numerically this property by fitting the time series for the ratio K(t) = 
E tr /E rot by means of the linear least-square fit K(r) = + a\T. We have performed a 
t-Student statistical test of the data to determine the 99% confidence intervals for a and 
ai— . The result of such a test shows that a\ is consistent with the value zero (see table II). 
As already established for mixture of inelastic hard spheres^ the ratio K deviates from the 
corresponding elastic value (K = 2), but attains a constant ratio during the homogeneous 
cooling regime. We observe from table II that the larger the inelasticity the larger the 
breakdown of energy equipartition, and the longer the dimer the larger the deviation. 

B. Velocity distribution functions 

We turn now attention to the velocity distribution functions. As discussed above one 
of the characterizing features of the HCS is the existence of a time independent scaling 
function Q(z). We have obtained numerically the related velocity distribution functions 
for the two translational components of the center of mass velocity and for the rotational 
velocity. These (normalized) distributions are displayed in Figs. El and compared with the 
Gaussian. We have not tried to measure departures from the Maxwellian, a task beyond 
the scope of the present work. Moreover, it is known from the literatur e) 1 ^ 30 that in the 
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case of inelastic hard-spheres these departures are small in the range 0.6 < a < 1 so that 
the Gaussian remains a valid approximation in the HCS of inelastic hard spheres. We guess 
from the present data that such a behavior remains true even in the case of the dimers. 

To summarize there is evidence that inelastic hard dimers display a normal behavior in 
the HCS even with respect to the nature of their velocity distribution functions which are 
nearly Gaussian. Thus the analogy with the hard-sphere system is complete. 
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FIG. 6: a) Normalized translational velocity distribution function for dimers of elongation a = d 
after different number of collisions per dimer (cpd). b) Normalized rotational velocity distribution 
functions for the same case. A gaussian curve (dashed line) is plotted for comparison. 



IV. CONCLUSION AND OUTLOOK 

We have investigated the dynamics of an undriven gas of inelastic hard dimers. The 
first achievement of this work has been to formulate properly the dynamics of pure inelastic 
dimer systems and that of a binary mixture of inelastic dimers and disks. Using such results 
we treated a pure gas of dimers with the ED algorithm. To validate the model and uncover 
some new physics we have then studied numerically the relaxation dynamics of such a system 
during the homogeneous cooling stage for various choices of the physical parameters. By 
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considering the total, translational and rotational kinetic energies of the dimers, we have 
found that their decay, following an initial short transient, agrees, as far as the system 
remains homogeneous, with the one observed in the case of a gas of IHS, that is Haff's 
cooling law. The different geometrical shape manifests itself in determining an identical 
pattern of cooling, although with a slower cooling rate which is related with the internal 
structure of the dimers. Moreover, we have found that the energy equipartition principle 
does not hold, although the ratio between rotational and translational energy is constant 
and almost equal to the equipartition value. 

The one-particle velocity distribution functions of the system remain always close to 
the Maxwellian with properly scaled parameters both in the case of the center of mass 
translational and rotational degrees of freedom. These observations allow us to conclude 
that a HCS exists also in the case of inelastic hard dimers. 

Preliminary work has also indicated interesting trends in the stages following the homoge- 
neous cooling regime, in particular we have observed that on a longer time scale the system 
breaks its homogeneity showing vortex structures in the velocity field (shearing instability) 
and clustering. Both phenomena have a counterpart in systems of spherically symmetric 
particles. However, a careful analysis of these results requires further work. 
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TABLE I: Dependence of the parameters D e and Vq*/vq on the coefficient of restitution and on 
the elongation of the dimer. For comparison the hard disk theoretical values are D e = 2 and as a 



numerical check 


we applied the 


same procedure as 


m the case of dimers and obtained D e 


= 1.998 


at area fraction 


7rnCT 2 /4 = 5.7 x 10" 4 . 










a 


D e 


V /^o 






0.99 


3.003 


1.07 






0.98 


3.013 


1.08 


a/d = 0.5 




0.95 


3.028 


1.08 






0.90 


3.079 


1.11 






0.85 


3.109 


1.13 






0.99 


2.987 


1.26 






0.98 


3.014 


1.25 


a/d = 1.0 




0.95 


3.038 


1.28 






0.90 


3.086 


1.32 






0.85 


3.041 


1.31 






0.99 


3.004 


1.45 






0.98 


3.031 


1.47 


a/d = 1.5 




0.95 


3.059 


1.50 






0.90 


3.050 


1.48 






0.85 


3.102 


1.50 
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TABLE II: Dependence of the ratio K between translational and rotational energy, on coefficient 
of restitution and elongation. The values Aao and Aai, concerning the 99% confidence intervals 
ao ± Aao an d a\ ± Aai, are also shown. 







a 


a 


Aa 


ai 


Aai 






1.0 


1.999 


1 x 10 -3 


-2.86 x 10" 6 


1.2 x 10~ 6 


a/a = 


0.5 


0.98 


2.014 


1 x 10~ 3 


—2.82 x 10 


3.5 x 10~ 6 








1 QQ4 


9 x in -3 


U.OO A ±U 


1 n v m~ 5 

X .VJ A 1U 






85 


9 056 


9 x 10 -3 


9 17x1 -4 


1 5 x 10~ 5 






1.0 


2.001 


1 x 10~ 3 


1.35 X 10~ 5 


2.1 x 10~ 6 


a/a = 


1.0 


0.98 


1.989 


1 x 10 


4.84 x 10 


3.1 x 10 






0.95 


2.046 


1 x 10~ 3 


-2.50 x 10~ 4 


8.6 x 10~ 6 






0.85 


2.064 


2 x 10~ 3 


-4.78 x 10~ 5 


1.7 x 10~ 5 






1.0 


1.987 


1 x 10 -3 


3.22 x 10~ 5 


2.0 x 10~ 6 


a/d = 


1.5 


0.98 


2.025 


1 x 10~ 3 


-5.38 x 10~ 5 


3.3 x 10~ 6 






0.95 


2.036 


2 x 10- 3 


-1.13 x 10~ 4 


9.2 x 10~ 6 






0.85 


2.075 


2 x 10- 3 


2.06 x 10~ 4 


1.6 x 10~ 5 
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